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Abstract 



A mechanism is proposed for initial stage of instability development that can induce the 
fragmentation of nuclear matter, arising as a result of collisions of non-relativistic heavy 
nuclei. Collision of heavy nuclei is simulated as a collision of two unbounded Fermi-liquid 
"drops". 

The instability origination in such a system is related to propagation of increasing oscil- 
lations in the nuclear matter. These oscillations can exist in a resting Fermi-liquid: modified 
Landau zero sound, modified spin and isospin waves, combination of these more simple waves. 
These instabilities are analogous to the beam instability in ordinary electron plasma. Behav- 
ior features of the obtained oscillation increase increments are provided. They can be used 
as indication for experimental confirmation of the proposed mechanism of fragmentation at 
nuclear collisions. Directions along which nuclear matter "jets" can be expected are specified. 

1 Introduction Berkley, Oak Ridge, USA; Hamburg, Germany; 

Orsay, Caen, France). This interest is caused by 
the prospect to confirm our notions about inner 
The fragmentation phenomenon, i. e. the simul- struc t U re of nuclei and the character of inter- 
taneous decay of the excited nucleus into lighter nucleon interactions and also by the possibility 
nuclei (fragments) and separate particles, has to obtain new fundamental knowledge about the 
been known for a rather long time. However origin of intranuclear and intranucleon forces, 
the scientific interest to this phenomenon have 

been rapidly growing lately due to the experi- Nuclear matter formed by the collision of 
ments on collisions of heavy fast nuclei, which heavy nuclei is an unstable object decaying 
were carried out in different scientific centers of within a very short time. The decay of arisen 
the world (CERN, Switzerland; Dubna, Russia; nuclear matter is accompanied by the forming 
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of new lighter nuclei and nuclear fragments (so- 
called nuclear fragmentation process). Nowa- 
days while describing the evolution of nuclear 
matter, formed after heavy nuclei collision, one 
usually considers two scenarios (see for example 
[H, [2]). According to the first one the state of 
statistical equilibrium can be reached in the nu- 
clear matter (the thermalization of the nuclear 
matter) permitting the thermodynamic descrip- 
tion of the fragmentation process by means of 
the phase transition theory methods. Accord- 
ing to the second scenario the resulting nuclear 
matter is essentially unstable matter, in which 
the processes allow only the dynamic descrip- 
tion. 

In the present paper we assume the latter 
scenario. Each colliding nucleus is a system of 
nucleons and it can be considered as a gener- 
alized Fermi-liquid, whose particles have spin 
and isospin (distinguishing protons and neu- 
trons) degrees of freedom. Let us note that the 
Fermi-liquid approach of Landau and Silin [3], H] 
is widely applied in theoretical nuclear physics 
in order to describe the properties of heavy nu- 
clei [5], [6] . In solving the specific problems of 
theoretical nuclear physics it is often assumed 
that heavy nuclei (with A > 100) can be con- 
sidered as the infinite drops of nucleonic Fermi 
liquid in cases when surface effects are unimpor- 
tant. Then such assumption enables us to sim- 
plify complicated mathematical treatment with- 
out violation of the final results. Note, that 
theoretical approach for finite Fermi liquids was 
elaborated by A.B.Migdal ( see [5] and refer- 
ences therein). 

In the present paper we also suppose the nu- 
clear matter, formed due to heavy nuclei colli- 
sions, to be the Fermi-liquid consisting of par- 
ticles with spin and isospin degrees of freedom. 
The inter-particle interaction in this kind of nu- 
clear liquid is described by such parameters as 
Landau amplitudes, by analogy with the normal 



Fermi-liquid theory [21 E]. So in order to con- 
sider the fragmentation formed after collision we 
come to the less complicated problem of interac- 
tion of two drops of nucleon infinite Fermi-liquid 
moving relatively to each other. Note that de- 
scription of collisions of heavy ions as two drops 
of infinite quark Fermi-liquids leading to forma- 
tion of quark-gluon plasma was analyzed in pa- 
per [7]. 

We have demonstrated that interactions of 
two nucleon Fermi-liquids in the system can 
originate a number of instabilities associated 
with increased oscillations of modes existing in 
the ordinary Fermi-liquid. Such waves include 
a modified Landau zero sound, modified spin 
waves, waves related to excitation of isospin de- 
grees of freedom, as well as more complex waves 
related to a combination of the mentioned sim- 
ple waves [H].The present instabilities are analo- 
gous to the known "beam" instability in the or- 
dinary electron plasma, occurring when charged 
particle beams pass through the plasma [8j [9]. 
Since the mentioned instabilities can evolve ex- 
actly after the mutual penetration of collid- 
ing Fermi-liquid drops and formation of united 
Fermi-liquid, we consider the formed system as 
an object which is far from the statistical equi- 
librium and requiring the dynamical description. 
The dynamics of such system in this work is 
described by the kinetic equation for the one- 
particle distribution function for nucleons in col- 
lisionless approximation, which is widely used in 
the theory of ordinary Fermi-liquid [3l 0] . Such 
approach for phenomena we consider is valid in 
the case when the characteristic time of insta- 
bility evolution is less in comparison with the 
characteristic time of relaxation in our system. 
In many cases such approach allows the ana- 
lytic description of the initial stage of instabil- 
ity development in formed nuclear matter with- 
out numerical calculations. Note that the use of 
mentioned kinetic equation restricts our consid- 
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eration to non-relativistic nucleoli energies. The 
modification of the kinetic equation for the de- 
scription of evolution of the nuclear matter with 
relativistic nucleons is given in JTOj [TT] . In these 
papers the development of filamentary instabil- 
ities in the nuclear matter at the stage of the 
mutual penetration of the two colliding nuclei 
was considered on the basis of the kinetic equa- 
tion for relativistic nucleons. 

Use of this approach allows us to describe 
the initial stage of instability development in 
nuclear matter, formed after the nuclear colli- 
sion and to obtain the dispersion equations for 
zero sound oscillations in such system (analogue 
of zero sound oscillations in ordinary Fermi- 
liquid). The solutions of equations for small 
Landau amplitudes are investigated in details 
in case of collision of nuclei in normal state as 
well as in case of collisions of heavy excited nu- 
clei. The expressions for the increments of zero- 
sound oscillations in nuclear matter are found. 
We have shown that those increments have char- 
acteristic logarithmic dependence on the excita- 
tion energy of nucleus when the nuclear matter 
is formed due to the collision of heavy excited 
nuclei. This characteristic dependence can be 
used for experimental proof or disproof of the 
mechanism of the initiation of instabilities in nu- 
clear matter which is suggested in present paper. 

We predict the directions (angular distribu- 
tion) in which it is possible to expect the matter 
jets, consisting of the nuclei fragments and sep- 
arate nucleons and originating as a result of the 
destruction of nuclei under collision. In this pa- 
per these propagation directions are related to 
the direction and value of the relative velocity 
of two colliding nuclei. The occurrence of such 
jets is caused in our opinion by the further evo- 
lution of instabilities, described in this paper. 
It is known that the experimental exploration 
of the angular distributions of outgoing matter 
jets is the important data source revealing the 



processes in formed dense blob of nuclear mat- 
ter. For this purposes the special experimental 
complexes are created, they allow to perform the 
mentioned measurements almost in the whole 
spatial angle in [2j [12]. The experimental de- 
tection of mentioned jets in predicted directions 
would by evidence of the realization of suggested 
mechanism of initiation of instabilities in collid- 
ing nuclei system. 

We suppose the confirmation of predictions 
given in present paper should be expected for 
example in experiments on nuclear collisions 
Gd+U or Xe+Sn (INDRA, [21 [12]) with incom- 
ing nucleus energy above 145 MeV per nucleon. 

2 Basic kinetic equations 
of two colliding Fermi- 
liquids 

A solution of the formulated problem on kinetics 
of two nucleon Fermi-liquids collisions is based 
on the assumption that the Fermi-liquid energy 
is a functional E (/) of one-particle distribution 
function /, (x, p, t) of quasiparticles (nucleons; 
here index V is used to denote a quasiparticle 
spin and isospin, x is a coordinate and p is a 
quasiparticle momentum, t is a time). We want 
to emphasize that in the Fermi-liquid theory the 
introduced energy functional plays the role anal- 
ogous to Hamiltonian in the microscopic theory. 

Evolution of one-particle distribution func- 
tion fi (x, p, t) in the collisionless approximation 
{ujT r ^> 1, to is a frequency, r r is a relaxation 
time) is governed by the kinetic equation 

dfi_ dei_ dft _ dtn dft = Q 
dt dp <9x <9x dp 

where quasiparticle energy e$ (x, p, /), deter- 
mining kinetics of the nucleon system, repre- 
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sents a variational derivative of the energy func- 
tional over the one-particle distribution function 



SEjf) 



A specific form of quasiparticle energy 
Si (x, p, /) as a functional of the distribution 
function in the Fermi-liquid theory is not known. 
Therefore, let us introduce functions of quasi- 
particle interaction F^ (generalized Landau 
amplitudes), which represent a linear reaction 
of the quasiparticle energy on small deviation 
Sfi (x, p, t) of distribution function /j (x, p, t): 



5ei (x, p,t) = 
= dr / c?x Fij ^ 

xdfj (x',p',t 



x-x ;p,p x 



(2) 



where dr = d 3 p/ (2iih) 3 . Functions Fij, repre- 
senting second order variational derivatives of 
the energy functional with respect to the one- 
particle distribution function, are main charac- 
teristics of the theory that can be experimen- 
tally determined (in this connection see, for ex- 
ample [T51 [T4"] ) . In view of short-range interac- 
tion forces between nucleons, from here on we 
will consider that the functions F^ have the fol- 
lowing form: 



x — x 



',P,p') = Fij (P'P ) 6 ( x 



(3) 



where quantities F^ (p,p ) (we will call them 
also Landau amplitudes) do not depend on co- 
ordinates. Note that assumption is valid 
only if we neglect the Coulomb forces induced 
by the proton charge. Linearizing equation 01) 
near a stationary quasi-equilibrium state (see 



below), described by distribution function f 0i , 
we can go over to Fourier components of devia- 
tions 5fi = fi — fo of the one-particle distribu- 
tion functions from their equilibrium values 

Sfi (p,w,k) = 



d 3 x / dtSfi (x, p, t) exp (iut — ikx) 



(4) 

As a result, using d2J) we get the following 
kinetic equation for 8f i (p, u;,k): 



J dr'^Fj (p,p)Wp,w,k) =0, 

(5) 

where amplitudes F^ (p, p') in accordance with 
(EJ), (j3J) are given by the formula 



Fa | p, p 



S 2 E(f) 



Sfi(p) $fj(p) 



fi=fi0 



and Em is the energy of noninteracting nucleons, 

p 2 

e = — (m is the nucleon mass). 
2m 

If we assume that interaction between nu- 
cleons is invariant with respect to the spin and 
isospin transformations, then a structure of Lan- 
dau amplitudes F for the nucleon liquid will be 
of the following form: 



F = J(»)/Wir(0) + j(<) (o-o-) 2rW+ 

+ (crcr) (tt) F^ (6) 

where F<®, F&, F®, are generalized Lan- 
dau amplitudes, I^ s \cr are unit matrix and 
Pauli matrices in the spin space, /W,r are a 
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unit matrix and Pauli matrices in the isospin 
space. Accordingly in equations ([SD there will be 
one of the Landau amplitudes from expression 
([6]), depending on which oscillations are consid- 
ered (oscillations of density, spin density, isospin 
density or spin-isospin density). Thereby, here- 
inafter we will deal with only one linearized ki- 
netic equation 

(^-kf^) 5f(P,w,k)+ 

J dr'F (p, p') Sf(jp, u, k) = 0, 

(7) 

where we will bear in mind F (p, p ) as one of 
the amplitudes from (jfJJ). 

Now let us obtain the expression for a quasi- 
equilibrium distribution function from kinetic 
equation ([7]). The equilibrium distribution func- 
tions of the resting drops of Fermi-liquid will be 
denoted by /g (p) and (p). Then in the 
case when the first drop is resting and the second 
one (impacting) is moving with velocity u, equi- 
librium distribution functions will be (p) 
and / (p - mu), respectively. This implies 
that the quasi-equilibrium distribution function 
of two colliding unbounded Fermi-liquids con- 
sisting of identical particles has the form 

/o (P) = «i/o (1) (P) + «2/ (2) (P - mu) . (8) 

Obviously that this function obeys the kinetic 
equation (CQ). 

We want to explain a meaning of coefficients 
ax and a 2 appearing in this expression. With 
this aim we will introduce auxiliary distribution 
functions g (p), gjp (p), g^ (p) normalized to 
unity (and because of that permitting a proba- 
bility interpretation) and related to each other 



analogously to relation (JSJ) : 

£o (p) = Qi9o (p) + Q2go ] (p - mu) . (9) 

From the normalization condition of these 
functions to unit it follows that values q\ and 
q 2 satisfy the equality 

qi + q 2 = 1. (10) 

From here we see that coefficients q\ and q 2 
must be interpreted as probabilities that an ar- 
bitrary selected particle belongs to either a sys- 
tem with distribution function g^ (p), or to a 
system with the distribution function g (p). 
However, in this case probabilities q% and q2 
must be specified by relative frequencies of the 
events of finding objects with distribution func- 
tions g^ (p) and gjp (p). If we are dealing with 
two colliding beams of nuclei, then for proba- 
bilities q\ and q 2 , we can obtain the following 
expressions (see ( TTUl) ) 

Ql = N x / (Nt + A 2 ) , q 2 = iV 2 / (Ai + iV 2 ) , 

(11) 

where N% and N 2 are densities of nuclei in the 
first and second beams respectively. 

As it was previously noted the above given 
argumentation implies the normalization of the 
distribution functions to unity. In the case un- 
der consideration this requirement is not fulfilled 
for functions f (p), / (1) (p), f^ 2) (p) because 

J dr fo (p) = J dr fo 1] (p) = n i> 

J drff (p) = n 2 , (12) 

where n is density of number of particles in the 
system consisting of two colliding Fermi-liquids; 
rii and n 2 are densities of number of particles 
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in each of these Fermi-liquids. But assuming in 
(ED that 



9o IP) 



/o(p) rf'(P) 



n 



fl (P) 

Tit 



™2 



and introducing denotations 

oc\ = nqx/n x , a 2 = nq 2 /n 2 , 
we can arrive at the expression (JH] 



(13) 



3 Zero-sound dispersion 
equation in a simple 
model 

In order to solve kinetic equation ([7j) we will use 
a simple model, where the Landau amplitudes 
do not depend on momenta p, p . We also as- 
sume that the colliding Fermi-liquids have zero 
temperature, so that the equilibrium distribu- 
tion functions of separate drops are given by the 
formulas 



<*/(p>w,k) (a; - kv) + 
dfo (P) 



+k 



dp 



F / dr6f(p,to,k) = 0, (15) 



where v = — = — , and in accordance with 
op m 

(El), (fH 



-a±v5 (eif — s) — a 2 (v — u) x 



dfo 
dp 



x 8\e 2F - 



m (v — uy 



(16) 



The most general solution of equation (TT5" 
looks as 



Sf (P,k,^) = 

= -F {to - kv + to}- 1 k ^°( p % (w, k) + 



<9p 



+<5A (p, k)8(u- vk) , 
where we denoted 

6<p (oo,k) = y c?r5/ (p,k,u;) 



(17) 



and 5 A (p, k) are arbitrary functions satisfying 
the condition (see (HI)) 



f (i) 
Jo 



p) = 0(e 1F -£),./o (2) (p)=0(e 2F -e) 
£o (p) , (14) 



g 

2m 



<L4* (p,k) = <L4(p, -k) 



(19) 



where 6(e) is theta function, e-yp, £if are Fermi 
energies of the colliding drops. It is evident from 
(fl4|) that in this case a difference between nu- 
cleus types can be found only in their differ- 
ent Fermi energies (note that for heavy nuclei 

£\F ~ ear). 

Within used assumptions linearized kinetic 
equation © can be written as: 



Formula (fl7l) permits to find a value of 5(p (to, k) 
in terms of the functions 5 A (p, k): 



5ip (to, k) = e 1 (to, k) 5 A (to, k) 
where (see flTQ]) ) 



(20) 



5A(to,k)= drSA(p,k)5(to -kv), (21) 
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8 A (uj, k) = 5 A (-to, -k) . 

Substituting further fl20l) into fTTTj) . we will 
have the following relation for 8f (p, k, co): 

<5/(p,k,^) = 

= 6 A (p, k) 5 (uj - kv) - F {to - kv + io}~ 

xg- 1 (u; ) k) ( L4(cu,k)k g/ ° (p) , (22) 

ap 

where e (a;, k) in (|2"UI) . f[2"2l is defined by the for- 
mula 



cj (k) + «7 (k) should be found from the equa- 
tion 

£ ~(u;o(k)+i7(k),k) = 0, (24) 

with decrement (increment) 7 (k) determined by 
the imaginary part of quantity i(uj, k). This is a 
reason why weakly damping or weakly increas- 
ing oscillations in the system 

Mk)|»| 7 (k)| (25) 
can exist only on condition that 



e (uj, k) = e* (-uj, -k) = e x (co, k) + ie 2 (uj, 

= 1 + Fk / dr — {co — kv + 10) 

J op 

(23) 

For the case of the charged Fermi-liquid 
i represents a complex dielectric permittivity 
(see, for example, [15]). As is known the pres- 
ence of an imaginary additonal component in 
e(u;,k) indicates damping or increase of wave 
amplitudes. The wave dispersion law co (k) = 



) \£i{u,k)\ > \i 2 (u,k)\. (26) 
Using the formula 

(z + io)~ l = P in5 (z) 

z 

(where P is a symbol of the principal value), 
and taking into account (fTB"]) . we will write func- 
tions £\ (uj, k), s~2 (k>, k) (see f l23|) ) in the follow- 
ing form: 



e~i (co, k) = l+oiiF < 1 



CO 



2kv 



In 



IF 



CO + kVip 



CO 



kv 



IF 



LO 



ku cos a 



2kv 



In 



2F 



CO 



ku cos a + kv 2 F 



CO 



ku cos a — kv2F 



e 2 (uj, k) = -F {ax- ( 

2 [ kvip 



CO 



kv 



IF 



V2F\UJ 

+ a 2 [ 

VlF, 



ku cos a 



kv- 



2F 



eli 



CO 



ku cos a 



kv- 



2F 



(27) 



Here T is a dimensionless Landau amplitude velocities V\p,V2F are given by relations 



mv 



£\F 



IF 



mv. 



£2F 



2F 



(29) 



Fm 2 V\p 



and angle a is an angle between directions of 

(28) 

wave vector k and velocity of impacting drop 
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u. In accordance with (1241) - (127)) a dispersion 
equation has the form e,\ (u (k) ,k) = 0, or 



1 + a x T { 1 - - In 



s + 1 



+ a 2 -Fx 
V 



x <! 1 ~ 2 ^ S ~~ ln 



s - 1 

77s - s + 1 



7]S - S - 1 



(30) 



and the decrement (increment) is given by 



7(*) 



dei (u, k) 
du 



-1 



e 2 (u,k), (31) 



7T 



e 2 (uj, k) = -T {a x sQ (1 - \s\) + 

+ a 2 ~ (r/s - s ) 6(1- \rjs - s Q |) 1 . 
77 J 

We use the following denotations in formulas 
©QD, (M) (see also ([28]), ([29])): 



_ u _v 1F _ u 

s = - , r] = , s = cos a (32) 

kvi F v 2 f v 2F 

(remember that for the heavy nuclei E\f ~ b 2 f, 
that is, 77 ~ 1) . 

dii (u;, k) 



The derivative 



in J3B), 



LU=LJQ 



according to (1271) . (13"0]) . acquires the form 



9£i (a;, fc) 



s 

u 



ln 



s + 1 




s - 1 


' s 2 - 1 



Solution of the dispersion 
equation for zero-sound 
oscillations 



A dispersion law of the zero-sound oscillations 
uj = u>o (k) is determined by the solution of equa- 
tion (|30)) . Note first of all that this equation 
(in the same way as the relation for j(k), see 
(13T)) . (1531 ) is invariant with respect to a simul- 



taneous substitution s 



s 



-s , (see 



(152)) ). This condition shows that propagation of 
the zero-sound wave is possible in opposite di- 
rections with the same increment or decrement. 
Consequently, let us consider s > for the sake 
of definiteness. Since equation (15U1) is quite com- 
plicated in its general form, we will make certain 
assumptions in order to obtain a solution. 



In 



7]S - S + 1 



rjs - s - 1 



+ ol 2 T 



(77s - s ) - 1 J 
(33) 



It is known that in the ordinary Fermi-liquid 
propagation of undamped zero-sound oscilla- 
tions at zero temperature is possible only under 
the condition JF > 0, s > 1 (in this connection 
see [ini UT] ) . Otherwise oscillations will quickly 
damp. Let us make an assumption that 



T > 0, s > 1 



(34) 



and demonstrate that in this case oscillations 
can increase or damp due to existence of the im- 
pacting drop. Naturally, one can easily see from 
(15T1) that possibility of propagation of weakly 
increased or weakly damped oscillations in the 
case, determined by inequalities (1541 . is related 
to the fulfilment of the condition 



[77s — So| < 1. 



(35) 
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If inequalities (13"1"1) . (|3"5|) are valid, quantity 
e 2 (t^Oj ^) = ^2 (s) has the form: 



7T 



£ 2 (s) = «2— (t/s - so)^ 7 . 
2?7 



(36) 



Let us study the solution of dispersion equa- 
tion (130]) with small Landau amplitudes J 7 , 
T <C 1 and at the conditions (134]) . fl35|) . Since in 
this case quantity s is close to unity, s ~ 1 + 5s, 
equation ( |30l) can be written as 



1 + a x T + —T - On - s ) x 
1] 2rj 



x In 



1 + (V ~ sp) 
1 - (V ~ s ) 



T^ ln 5? 



whereof 



5s 



1 + (V ~ s ) 1 ^ ( "~ So) 



1 - (?7 - s ) 

«2 



xexp -2 1 + 



OtlT) 



exp 



(37) 



where 



s = 1 + 5s. 

Taking into account (1331) . for conditions 
f l35|) we have 



dii\ IT . 

OUJ I , , Un 25s 



Further using formulas (|3~TT) , (|36|) , (!37|) , we have 
arrived at the following relation for 7 (A;): 



— — PS! 27T (77 — S ) X 

^0 

1 + (V ~ so) 



1 - (V - s ) 
2 



L (»7-so) 



X 



x exp 



1 + ^^ + 



(3J 



1 77 — s 1 < 1 , uj = kvi F , .F<Cl. 



According to (j4]), ([21]) the condition 
7 (A;) > 0, or < (77 — s ) < 1, corresponds to 
the oscillation damping, whereas the condition 
7 (A;) < 0, or — 1 < (77 — s ) < corresponds to 
the oscillation increase, meaning instability de- 
velopment. The latter inequality can be written 
(in view of (132]) ) as 77 < So < V + 1> ° r taking 
into account rj = \/e\fI£2F 



V\f < UCOSa < V\f + V2F, 



(39) 



where a is an angle between vectors u and k. 

It is evident from (1381) that when approach- 
ing to the right boundary of interval (|39|) . that 
is if 



s 



77 + 1 , s = l + 5s 



(40) 



the increment 7 (k) is unrestrictedly increasing. 
However, it is evident that for the case PUI) we 
need to solve more correctly dispersion equation 
( 1301) . which in this approximation and keeping 
in mind s — > 1 has the form: 



«2 r 1 sr- 1 5s a 2 . r/5s 
77 2 2 zt) 2 

The solution of this equation is given by the re- 
lation: 
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a 2 



5s 

x exp 



2r) Vi a i + «2 x 
-2rj 



1 + a\T 



V J 



(77CK1 + 02)^ 

Remembering that according to (1331) un- 
der conditions (1ST]) . (HOI) the relation for 

So; J 



is 



LJ=UJq 

dei 



1 r/ai + a 2 T 
26s 



_=u,o 77 

and considering fl31~]) . fl36|) . (j4]), we can obtain 
the following formula for increment 7 (/c): 



Fermi-liquid drops have zero temperature (see 
(1141) ). However, when two equilibrium Fermi- 
(41) liquid drops collide and the temperature of each 
one is non-zero, a scenario of weak instability 
development can have significant peculiarities. 

In order to solve this problem we need to 
modify relations f[5U|) . (|31|) . determining fre- 
quency and increments (decrements) of zero- 
sound oscillations in conformity with the case 
of the colliding equilibrium Fermi-liquid drops 
with the non-zero temperature. In this case the 
equilibrium distribution functions of the drops 
are as follows (see (jSJ) - ffT3l) ): 



7 (*0 
x exp 



a: 2 



-2tt- 



«2 



r]ai + a 2 
-2rj 



7/ 



ijai + a 2 



x 



(77a 1 + a 2 )T 



1 + ai.F + 



(42) 



.F<1, u &kv 1F , s«l, s < 77 + I, 

where quantities sq, tj, s are defined as be- 
fore by relations (|32|) . Taking into account that 
(rjai/rjai + a 2 ) < 1, comparing relations (1551) 
and (42) one can easily conclude that when 
T <C 1, the maximum value of the increment is 
reached for the upper limit of the interval fl39|) . 
i.e., when 



V^E cos a < v/^If + v^2F' (43) 

where £" is the kinetic energy per one nucleon 
in the impacting drop. 



5 Temperature impact 

So far we have discussed the problem of prop- 
agation of weakly increasing or damping zero- 
sound oscillations, assuming that the colliding 



tf>(p) 



/o 2) (P - mu ) 



exp 



2m 



— E\ F 



+ 1 



exp 



(p — mu) 

2m 



^2F 




(44) 



where T 1; T 2 are temperatures of the resting 
and impacting drops, respectively. The given 
below relations should hold true for the drops 



(Ti/e lF ) < 1, (T 2 /e 2F ) < 1. 



(45) 



The relations determine possibility to apply 
the Fermi-liquid description to systems of strong 
interacting fermions (in the present case - nu- 
cleons). Taking into account formulas (jHJ) and 
inequalities (45), one can make low temperature 
expansion in expression (123|) . that can result in 
the following equation for determination of fre- 
quency ujq = skvip of the zero-sound oscilla- 
tions in the main approximation over parame- 
ters (Ti/ei F ), (T 2 /e 2F ) (see [H]): 

£i(s) =£io(s) +e 1T (s) =0, 



10 



- Ml l.sj 1 • a x T{ 1 - - In 



s + 1 



l + a 2 ^-T jl - ~ (rjs - s ) In 



rjs - s + 1 



7]S - s 



e 1T (s 
— a\T 



P 



zdz 



2 + (s 2 -l) 



e + l ( S 2 _ i)^ ( £lF / Tl 



a 2 .F 



P 



zdz 



2 + [(r/s 



-i] 



where P means as before the symbol of the prin- 
cipal value. In the main approximation the ex- 
pression for quantity £ 2 (s), defining in accor- 
dance with ff271) . (l3Tj) a weak increase or damp- 
ing of the zero-sound oscillations, can be put as: 



«2 - S ) 



V 



fo ( E -f 2 [{vs - s ) 2 - 1] 



(47) 



where 



f (e/T) = {e s / T +iy 



(4£ 



It is easy to see that at Ti — > 0, T 2 — * quan- 
tity Eyr (s) turns to zero. Hence, equation (1461) 
transforms into dispersion equation (1301) . Ex- 
pression ()47l) . with account of ()48l) . transforms 
to formula ( 13TT) for quantity £2(5), determin- 
ing instabilities of oscillations when the drops 
have zero temperature. As we can see from (|46|) . 
f )47|) . the temperature effects can greatly impact 
development of the instabilities associated with 
propagation of the zero-sound oscillations in the 
system, when s ~ 1, i.e. J 7 < 1. Keeping in 
mind this fact, let us analyze the possibility of 
weak increase of the zero-sound oscillations in 



+ 1 [(rjs - s ) 2 - 1] (e 2F /T 2 ) 2 - 

(46) 



the system when the resting Fermi-liquid drop is 
found at zero temperature, T\=0. Here we will 
also assume that the following relations (com- 
pare with (|M]), (Hi, (H2), (PD, (45)) are valid 



.F<1, s>l, s < 77 + I, 0WT 2 )>1, 



5s {e 2F /T 2 ) < 1, Ss = s - 1 « 1. (49) 

Under these conditions dispersion equation 
meeting these conditions, in the main ap- 
proximation over Ss will get the form: 



Ei (s)«l + at! T + a 2 ^ (1 - Jo) - 



a 2 T , 4e 2F uiF , Ss 

In 1 In — 

V T 2 2 2 



0. 



(50) 



In order to obtain this equation we used 
the asymptotic estimate of the integral, when 
1*1 < 1 



P 



zdz 



(c + i) (t2-^|7^ /o+ 2 ln|t| ~ Jlt ' 

(51) 
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1 

4 
i 



dz 



z 



i 



1 f dz r z z 

H — / — < tanh 

2J z \ 2 2 

o 



0,07, 



h 



dz 



z 3 (e z + 1) 



dz 

z 



. z z 
I tanh 2"2 



we can arrive at the following expression for in- 
crement 7^ (s), considering (1511) . 



x exp 



a 2 

7T— — te^X 

1 r i + ai^+^(l-/ ) 



»7 



0,11. 



(55) 



which is true at the conditions (1491) . (1531) . In 
order to obtain formula (EH), we have used ex- 



The integrals having exactly the same form de- 
scribe the temperature impact on dispersion Passion flE), and also have taken into account 

equation gg]). The solution of equation Q5DD is: that ( see ® ) 



= 2 ( 4 f|) ^ x exp {_^} , (52) A (f " 0) - o (i - s 2 H o, . > 1, 



Qi = l + ai^+ — (l-/o)- 

In accordance with (T491) this solution must sat- 
isfy relation 5s (^f/^) <C 1, whence, remem- 
bering that (S2F/T2) ^> 1, it is easy to deter- 
mine a "temperature" range, when solution (l52l 
exists for 5s: 



So 77 1 

K — <C exp 

J-2 



277Q! 



(770*1 + a 2 ) ^ 



(53) 



Noting further that under (150]) the formula 
holds 



dei 



1 

cuq 25s ' 



lu = kV\F 

and using the fact that when T\ — > and s > 1 



fo 



£ 2F , n2 1 

— (77s - S ) - 1 

J-2 



/o(0) 



1 

2' 



so < ^ + 1- 

Now let us consider a case when the im- 
pacting Fermi-liquid drop has zero temperature, 
T 2 =0, and the resting drop has temperature T 1 
We assume that conditions (compare with ( 149]) ) 
are valid: 



JF<1, s>l, s < 77 + I, (e 1F /T 1 )>l, 



4 77 



(54) 



5s (Siir/Ti) < 1, fc = s-l«l. (56) 
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In this case the dispersion equation in the main 
approximation over 5s considering asymptotic 
estimate floTl) has the form: 



the expression for e 2 , determining in accordance 
with (13 ip the damping or increase of the zero- 
sound oscillations, can be put in a form 



ii(s) « 1 + 01^(1- J ) + 



ol 2 T 



71 



e 2 = —T{ria\ - 2o 2 } . 
4?7 



(60) 



Bearing in mind that quantity {dei/du} 



UJ = UJQ 



In -^7— + ^ n = 0- (57) in compliance with (1571) is given by the formula 



A solution of this equation is given by 



63 = \ ( 4 ^f) ^ X ^ 



Q 2 = l + ai F{l- J ) + 



2r?Q 2 

2 JF 



(5f 



o 2 JF 
77 



and here, we should remember that in accor- 
dance with ([55]) the relation 5s{e\p/T-\) <C 1 
must be valid, and a "temperature" condition 
of existence of such solution is determined by 
the inequality (compare with (153]) ): 



1 „ E-if ( 2i]Q 2 

1« — «exp ■ — : 

Ti ^ (770:1 + o 2 ) T 



(59) 



dei 
duj 



1 a 2 T 1 



(61) 



=wo ^0 2r/ 5s 

where 5s is determined by fl58|) . and substituting 
(EOD, (JMD into (EU), we have 



Tfc 
^0 

x exp 



7T- 



77«i 



2o 2 / 4 fi£ 
t?o 2 V Ti 



o 2 jF 



1 + 01^(1-/0) + 



ol 2 T 

V J 

(62) 



UJQ 



kv 



IF- 



Noting further that for s > 1, So < r] + 1 and 
T 2 = the following formulas are valid 



/o(f (3 2 -l))-/o(2^)-/o(0) = l 



/o [fa - *o) 2 - 1]) ^ 6 (1 - (r/s - s ) 2 ) 



As is easy to see, the formula (|62|) in contrast 
to (|55j) . which is correct when Ti = 0, T 2 7^ 0, 
determines either a coefficient of damping (when 
r/01 > 2a 2 ), or increase (when 770:1 < 2o 2 ) of the 
zero-sound oscillations. Such non-invariance of 
fl55l) . (162]) relative to interchanging of the im- 
pacting and resting drops with a correspond- 
ing change of their characteristics e\p <-> £ 2 f, 
T\ <-> T 2 (in other words the absence of the 
Galilean invariance in the system of the collided 
drops) is not in any sense paradoxical. This is 
because within the developed Fermi-liquid de- 
scription model the interaction of quasiparticles 
is not invariant with respect to the Galilean 
transformations (see ([2]), ([3])). 

Therefore, the relation 7701 < 2o 2 that can 
be written as 



fo 



£2F 

7^ 



[(77s 



•soj 



-1) 



T 2 ^0 



iE\p < 2o 2 -i/£ 2 F, 



(63) 
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along with (j56j) . (j59p determines existence con- and coefficient 1\ is determined by (15 ip . 
ditions of weakly increasing oscillations with in- The solution of equation (I6"6"j) . since 
crement 7& in the system {A/ B) <C 1, will be unambiguous and positive 

when A > 1: 



mi 



7^ r\a x - 2a 2 ( .e 1F 
rja 
2?/ 



»7«a V T J < y s = _ + ^__ + __ A>1. (68) 



x exp 



a 2 J = I r] 

(64) 



'Focusing on the explicit form of coefficients A 
and B, it is easy to prove that the conditions 



In order to finish studying the temperature ^ . ^ 

impact on the instability development in the S ^ 6lF ' ^ ' s ^ 2i? / 2 ' 

system of two collided Fermi-liquid drops, we are automatica ii y fulfilled. 

want to obtain an expression for the increment p rocee ding from fl3T]), © and noting that 

when both the resting and impacting drops have in accordance wit h (@ZD, dS5H, dSSD the following 

comparable (Tf ~ T 2 ) nonzero temperatures, f ormu l as hold 



moreover, 



So ~ J 1 ' iv a l — «2i 

^«l, S >l, So <r ? + l, 477 
(£if/Ti) > 1, (£ 2 f/T 2 ) > 1, 



(A-l), 



. , /rTn . _ . , ;rTn . we will attain the following expression for coef- 

5s(e 1F / Tl )«l, 5s(e 2F /T 2 )«l, ^ 

fe=s-l<l. (65) 

Dispersion equation fT46l) in the main approx- ^ k n rja\ — a 2 

imation, with taking into account conditions Uq ~ Arj^ A 2 + AB (A — 1) ^' W ° ~~ t ' 1 ' F ' 



( 161)]) . has the form: 



(69) 



D , 2 /ir , 1 4 n /**\ It is clear that when 77a i < a 2 , or 

£1 = Bos — Ads + 1 — A = 0, (66) 



where 



jF f e 2 F 1 quantity 7^ represents an instability increment 

= "2" 1 ai ^ n ° 2 ^ n j 1 J ' of the zero-sound oscillations. 

Expression (169]) takes a simple form 

fi-wJ^f^y + ^feV) (67) Tl^JLj: ^ ~ 012 (71) 
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or 

— w ^-T (rjax - a 2 ) , = kv 1F 
when A ~ 1, that is (see (1S"7|) ) 

ai In (e 1F /Ti) + a 2 In (e 2F /T 2 ) ps -. (72) 



For a more general case, i.e. when A ^> 1 
and A i? (see (E7j)), formula (1691) also gets 
simpler, but it is still more complicated com- 
pared to (1711) : 



o; 



7T (j7«i 



«2 



3277^ 



[ai In {eip/Tx) + a 2 In (e^/^a)] 



-1/2 



ai 



+ r?a 2 



f)' 



Expressions for increments (1711) . (1731) look 
even simpler, if we consider collision of identical 
nuclei having the same temperatures (see be- 
low formulas (1751) , (1761) ) . Comparison of expres- 
sions ([33), (@2D with (USD, (|S2p, flTJj), JED, with 

taking into account relations (16"3"1) . (j70j) allows 
to come to a conclusion that for the colliding 
Fermi-liquid drops with nonzero temperatures 
the instability increments are higher compared 
to the ones for the colliding drops with zero tem- 
peratures. In particular, for the collision of the 
Fermi-liquid drops having comparable tempera- 
tures Ti ~T 2 ((Ti/eu) < 1, (T 2 /e 2F ) < 1), 
when condition (172]) is fulfilled, the increase 
increment of the zero-sound oscillations has a 
power like character of dependence over the Lan- 
dau amplitude T <C 1, whereas for other consid- 
ered here cases (except the one leading to (JT3"]) ). 
the degree of the increment smallness is given by 
exponential multipliers of the type exp (-X/J 7 ), 
T <C 1, where A is a certain constant. If we 
recall that in our consideration collision of two 
fast (but not relativistic) nuclei is simulated by 
collision of the two Fermi-liquid drops, then the 
provided comment means that the instabilities, 
related to increase of the zero-sound oscillations 
in the system of the two collided excited nuclei 



(73) 



(i.e. nuclei with nonzero temperatures), must 
develop much more intensively compared to the 
instabilities in the system of the collided unex- 
cited nuclei. 



6 Discussion of results 

Up to now we have studied the case of small 
Landau amplitudes f < 1. As it was proved 
in this case the dispersion equations of the zero- 
sound oscillations allow for an analytical solu- 
tion in the perturbation theory over small T . 
Analytical solutions of dispersion equations can 
be obtained also for the case of large Landau 
amplitudes, T ^> 1. However, if we consider 
these solutions from the viewpoint of the prob- 
lem about instability development in the sys- 
tem of collided non-relativistic nuclei, it would 
be easy to prove that a solution of the disper- 
sion equations in the perturbation theory over 
a small parameter 1/jF will lead us beyond the 
non-relativistic approximation, which was one 
of the main assumptions in the present work. 

Naturally, it is easy to prove that the solu- 
tions of dispersion equation (130|) are within the 
range of large s, s ^> 1 (note at once that when 
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s 3> 1, it is possible to neglect temperature de- 
pendent components in equation (T4"6"j) ). But in 
accordance with fpTTl) . f[33l) . f|471) for existence 
of instabilities the following condition must be 
satisfied so > or 



(1711) takes the form: 



yejep cos a > s ^> 1 , 



(74) 



where e is kinetic energy per one nucleon in the 
impacting drop. In inequality ( 1741) we took into 
account that for heavy nuclei eip ~ £ 2 f = £f, 
meaning, jj « 1. Remembering that for the 
heavy nuclei e F ~ 36 MeV, and the non- 
relativistic approximation requires the nucleon 
kinetic energy to be small compared to the rest 
energy mc 2 ~ 1 GeV, (e/mc 2 <C 1), it is easy to 



see that relation (I74|) can not be satisfied under 
the non-relativistic approximation. Therefore, 
in this work we do not consider the case of large 
Landau amplitudes. 

In this article we have restricted ourselves 
by consideration of positive Landau amplitudes. 
When Landau amplitudes are negative, the dis- 
persion equation of the zero-sound oscillations 
can be solved only by numerical methods. This 
is because in this real and imaginary parts 

of the frequency are comparable with each other 
and the analytical methods, based on the per- 
turbation theory, become inapplicable. The nu- 
merical solution of the dispersion equations for 
the negative Landau amplitudes is not included 
into the present work, which has as its object to 
demonstrate a principal possibility of instability 
development in the system of the colliding heavy 
nuclei, induced by propagation and increase of 
the zero-sound oscillations. 

Formulas fl7T|) . (1731) could be used SIS db 
proof of possible development of just such in- 
stabilities in the system of the heavy collid- 
ing nuclei. When two identical nuclei collide 
(sip = £2F — £f ,V = 1) T\ « T 2 = T), formula 



7^ ^ 7r ai-a 2 1 
uj 2 oc\ + a 2 In (e F /T) 



ujq = kv F . (75) 



In this case the expression (1731) is signifi- 
cantly simplified too: 



Ik 
x 



7T 



«1 — «2 



(e F /TY 



32 Ii \a l + «2 
{(1/2) T (a a + a 2 ) In (e F /T) - iy 1/2 , 

(76) 

As it is we known temperature T can be re- 
lated to the excitation energy of nucleus U by 
an approximate formula (in this connection see, 
for example, [HI M, ED [22] ) 



T 



all, 



(77) 



where a is a certain constant, which should be 
determined from experimental data (see in par- 
ticular [T9], [201 EI])- F° r example, in accor- 
dance with [19] for a nucleus with mass number 
A = 115, a « l/8MeV;A = 181, a w 1/10 MeV. 
Note that within the framework of the Fermi- 
gas model of nucleus quantity a can be found 
from formula 

Ae F 

Following this formula for e F ~ 36 MeV and 
A = 115, A = 181, we have a w 0, 124 MeV and 
a ~ 0, 0795 MeV, respectively. 

If we consider that the instabilities, associ- 
ated with the increase of amplitudes of the zero- 
sound oscillations, are responsible for fragmen- 
tation, then specific time of fragmentation ry 
must be of the order of the value reciprocal to 
increment 7^ 



Tf ~ l/7 fc . 



(7f 
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Thus, if we could experimentally mearsure the 
fragmentation time of nuclei Tf, then the specific 
dependence of the time from excitation energy 
U would be governed by the formula 



Tf 



Cln 



£ F 



c 



2 («! 

UJQ — 



«2j 



'aU J ' n (ax- a 2 

when relation (1721) holds, and by the formula 



(79) 



Tf 



c 



JaU 



T («i + a 2 ) In 



AZ7 



- 2 
(80) 



when condition A> 1, leading to (173]) is valid. 
According to (fT5l) - (178l) this would show that the 
proposed mechanism of the instability develop- 
ment in the system of the colliding heavy nuclei 
is working. Besides, if such instability develop- 
ment mechanism is realized, it can result in nu- 
cleus fragmentation and consequently, according 
to ([72D, ([79]), (HDD we could estimate both val- 
ues of frequencies of oscillations of developed in- 
stabilities and values of the Landau amplitudes 
for the nuclear matter formed as a result of the 
collision. It is should be noted that the Lan- 
dau amplitudes are parameters of the Fermi- 
liquid approach of the present work and, gen- 
erally speaking, they don't necessarily coincide 
with the Landau amplitudes for certain heavy 
nuclei, values of which are contained in some 
published papers (see, for example [TBI [Lf] ). 

However, if we make an assumption that val- 
ues of the Landau amplitudes in the nuclear 
matter, formed as a result of the collision of two 
fast heavy nuclei, are close to the correspond- 
ing values of amplitudes of certain nuclei, then 



we can estimate values of increments for spe- 
cific cases and provide suggestions about range 
of use for the results of the present work. Work 
[H] gives grounds to make a conclusion that 
the value of Landau amplitude Fq°\ determin- 
ing oscillations of the nucleus density (see §6§), 
is within the range of -0,25 and 0,5. These val- 
ues of the Landau amplitudes from the positive 
part of this range, when E\f ~ £2F = 36 MeV, 
are fully within the frame of the assumptions 
made in present work. 

The values of amplitudes F , Fq^ (see ([6])), 
determining oscillations of the spin and isospin 
densities in the nucleus (see [133 E]) are also in 
compliance with the assumptions of our work: 



, „ ~ 0, 27 -=- 0, 35, F (i) « 0, 59 -=- 0, 72. 

As for the negative values of amplitude F 
and values of amplitude Fq S ^ (see [13], UJ] ) 



F (si) « 1,26-=- 1,36, 

determining zero-sound oscillations of the spin- 
isospin density, these cases as we have already 
said require numerical calculations. This is the 
reason why they are not included into this work. 

Finally, let us study the question about pe- 
culiarities in behavior of a non-equilibrium dis- 
tribution function of two colliding nucleon drops 
in the momentum space. This will allow to de- 
termine a direction of nucleus fragment release 
result of instability development in such a 
system with respect to a direction of velocity of 
the moving drop. 

It follows from (]17l) . that deviation of dis- 
tribution function 5f (p, k, u) from the equilib- 
rium distribution function is determined by the 
equation 
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a i 



dti 1] (P) 
dp 



kv + io 



Q'2 



d/p 2) (P - mu ) 
dp 



10 



kv 



zo 



X 



X 



dr'S~f ( p, k, tu ) +SA (p, k) 5 (o> - vk) 



51) 



<9/o (1) (p) C(P-™ 



where derivatives 



(2) 



dp 



arc 



<9p 

given by (|T6|) . According to (jlj), (IHT!) the non- 
equilibrium distribution function can have max- 
imum, given the following conditions are met 

(^-kv)| v=ViF = 0, 
(^-ku-k(v-u))| |v _ u| _ 2F =0, 
v = p/m. 

Since to ~ c^o = skvip, these conditions can be 
written as 

s — cos 9 = 0, r/s — s — cos /3 = 0, (82) 

where 9 is an angle between directions of vec- 
tors k and v, f3 is an angle between directions 
of vector k and vector v — u. The first of these 
conditions is not fulfilled because s > 1. As 
we have previously mentioned a damping or in- 
crease of the zero-sound oscillations can take 
place only when \rjs — Sq\ < 1. This means 
that the second condition in (182]) can be ful- 
filled. The maximum value of increment 7& is 
reached when r/s — s ~ —1 (see fHU]) ). Under 
condition E\f ~ £2F = Sf ij] = 1) and s > 1 it 
corresponds to the condition 



u 

s = — cos a 

v F 



(83) 



(a is as before an angle between directions of 
vectors u and k) and, hence 



Relation r/s — s ~ —1 is accordance with 
means that vectors k and v — u are anti- 



parallel 



u 



or, bearing in mind that |v — u| = Vf, 



v — U = —Vp — . 

k 

Multiplying scalarwise this equation by vec- 
tor v — u, with taking into account ( l83il . we can 
obtain the following relation for cosine of angle 
9q between directions of vectors v and u 



cos 9n 



u 2 - 2vl 



vu 



Hence, remembering that the module of vec- 
tor v is within the region of values u — vf and 
u + vf, and also that by (JHl 



u 



2v 2 F 



u[u + vf) 
we have 



< 1, 



u 2 - 2v 2 F 
u(u — vf) 



> 1, 



-2vl „ „ 

-^r < COS# < 1- 



U (u + Vp) 

Taking into consideration the inequality ( 1841 
the latter equation can be written in the follow- 
ing form 



u > 2vf- 



34) 



< cos#n < 1- 



15) 
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Thus, according to the results of this paper 
in the case when the inequality (1831) is true the 
jets of nuclear matter at the collision of heavy 
nuclei should to be expected along the directions 
defined by the condition (185|) . In the case when 
the inequality (1831) is breaking, the angular dis- 
tribution of the outgoing matter should be close 
to an isotropic one. 

Let us come back to the problem collision- 
less approximation for the description of dynam- 
ics of nuclear matter formed by the collision of 
heavy nuclei. The main condition of its appli- 
cability ([1]) from the point of view of the results 
of this paper is the inequality loqt t ^> 1, where 
ujq = skvp and r r is a relaxation time. It is 
necessary also that the characteristic times of 
the instabilities' development tj ~ I/7 (where 
7 are the increments found in this paper) should 
be smaller or have the same magnitude as the 
relaxation time - < r r . The latter equation can 
be written down like — UnT r > 1. Thus, the 
equations 

U T r » 1, UJ T r > 1. 

^0 >> 7 > l/r r , t « I/7 < r r . 

are the applicability conditions of the collision- 
less approximation for the describing the initial 
stage of development of the instabilities in the 
nuclear matter, formed by the heavy nuclei col- 
lisions. Let us remark that the kinetic equa- 
tion without taking into account collisions be- 
tween particles had been used by some authors 
for the description of the dynamics of the mat- 
ter formed by the heavy nuclei collisions (see in 
this case [TJ EJ, [23j , and references therein). As a 
rule this equation was written in the mean field 
approximation. For example the spinodal fission 
in the expanding nucleon Fermi-liquid being was 



investigated in [23] on the basis of such kinetic 
equation. However, as it is known, the aver- 
age field approach needs considerable numeric 
calculations. The usage of the Fermi-liquid ap- 
proach with phenomenological parameters of in- 
teraction (Landau amplitudes) for the describ- 
ing the dynamics of the nuclear matter allows 
in many cases to use analytic treatment with- 
out the help of numerical methods. 

As it was repeatedly noted in this paper the 
application of Fermi-liquid approach is justified 
while describing the properties of heavy nuclei. 
For this reason the confirmations of the mecha- 
nism of origin of instabilities in the nuclear mat- 
ter suggested in this paper, one can expect in 
experiments at the collisions of heavy fast (but 
non-relativistic) nuclei. These might be the re- 
actions Gd+U or Xe+Sn (INDRA, E2]) at 
the energies of the incoming nucleus more than 
145 MeV per nucleon. 

The work is performed under STCU (Grant 
N 1480) funding support. 
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